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Abstract: In the study of the reaction —>• KsKi at the VEPP-2M collider with the 

SND detector the nuclear interaction length of meson in Nal(Tl) has been measured. Its value 
is found to be 30-50 cm in the Ki momentum range 0.11-0.48 GeV/c. The results are compared 
with the values used in the simulation programs GEANT4 and UNIMOD. 
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1. Introduction 

In the analysis of experimental data in the high-energy physics the Monte-Carlo (MC) simulation 
is widely used. There are known simulation packages, such as UNIMOD |[I]], GEANT3 [§], and 
GEANT4 ijUl, which successfully describe the physical processes of particles passage through the 
detector material. Nevertheless, in some cases, the accuracy of the theoretical models and/or the 
experimental data used in simulation may be not sufficient. Therefore, it is worthwhile to verify 
the simulation using new experimental data. 

In annihilation into hadrons at the center-of-mass (c.m.) energies above 1 GeV the 
neutral kaons are intensively produced, e.g. in the process 

e+e-^KsKf ( 1 . 1 ) 

Whereas Ks meson decays rapidly via Ks — )• Tin channel, the neutral Kl meson, due to its much 
greater lifetime, passes through detector material and is absorbed due to the nuclear interaction. 
Since kaons in the process ( [01 ) are producing in pairs, the Ki meson can be tagged by the recoiled 
Ks meson. This greatly facilitates the study of Kl interaction with the detector material. 

In this paper we use kaons from the reaction (jOJ) to measure the Ki nuclear interaction length 
in Nal(Tl) crystals in the momentum range from 0.11 to 0.48 GeV/c. 

2. Experiment 

The experiment on the study the —)• KsK^ process was carried out at the VEPP-2M [Q] 
collider with the SND detector |Q] in the c.m. energy range ^/s = 2E =0.98-1.38 GeV. The SND 
is a general purpose non-magnetic detector for low energy -colliders. In the center of the 
detector around the collider beam pipe a tracking system consisting of two drift chambers is in¬ 
stalled. The main part of the SND is a three-layer spherical electromagnetic calorimeter based on 
1640 Nal(Tl) crystals with the total weight 3.6 ton. The calorimeter covers the polar angle range 
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Figure 1. The distribution of y} of the kinematic fit to the Ks —?> \ j hypothesis for candidates 

in data (points with error bars) and —>• KsKt{Y) simulation (histogram) at = 1.04— 1.10 GeV. 


18° < G < 162°. The crystal angular dimension is A6 = A0 = 9°. The radial positions of the 
calorimeter layers are shown in Fig. The total calorimeter thickness for particles originating 
from the interaction region is 34.7 cm (13.4Xo), where Xq is the radiation length equal to 2.6 cm 
for Nal(Tl). Pairs of crystals of the two inner layers with thickness of 2.9 and 4.8Xo are sealed in 
a common thin (0.1 mm) aluminum container. In the gap between the second and third calorime¬ 
ter layers an aluminum supporting hemisphere is located. The scintillation light signals from the 
crystals are detected by vacuum phototriodes. Outside the calorimeter an iron absorber and a muon 
detector are placed. 

In this paper the data collected in the energy interval ^/s =1.04-1.38 GeV with an integrated 
luminosity of 9.2 pb^' are analyzed. These data were used previously to measure the —>• KsKl 

cross section [^. In addition, the data recorded at the maximum of the (^(1020) resonance with an 
integrated luminosity of 0.5 pb^' are used. 

The Monte-Carlo event generator for the —>• KgKi reaction used in this analysis includes 

radiative corrections [^. In particular, an extra photon emitted by initial electrons is generated 
with the angular distribution modelled according to Ref. |Q]. The —)■ KsKl Bom cross section 

needed to calculate the radiative corrections is taken from Ref. [Q]. 

The response of the SND detector is simulated using the UNIMOD package ||i|]. This package 
developed in BINP (Novosibirsk) in eighteens is functionally similar to the GEANT3 [§] simula¬ 
tion program. Special attention in UNIMOD is given to simulation of pion and kaon low-energy 
nuclear interactions. Nuclear cross sections are calculated with the SCATTER program [^, while 
the nuclear reactions are simulated using the NUCRIN model []E|. The simulation takes into ac¬ 
count the variation of the detector and accelerator conditions (dead electronic channels, size of the 
interaction region, etc.) and beam-induced background photons and charged particles overlapping 
events of interest. 


3. Selection of e^e —)■ KsK^ events 

The Ks mesons from the process KsK^ decay inside the SND tracking system. The decay 
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Figure 2. The Mrec spectrum for selected e+e — ?> KsKl candidates in four energy intervals: (a) -/s = 
1.04 - 1.05 GeV, (b) = 1.06 - 1.09 GeV, (c) = 1.10 - 1.20 GeV, and (d) y^ = 1.20 - 1.38 GeV. The 

shaded histogram shows the (on^ background. The vertical line indicates Kl mass. 

mode —)• —)• 47 is chosen for Ks reconstruction. Therefore we select events containing at 

least four photons and no charged particles. To eliminate the cosmic ray background, the veto from 
muon detector and the requirement that the fired calorimeter crystals do not lie along a straight line 
are used. 

We combine four photons in an event to create a Ks candidate. Photons included into the Ks 
candidate must be in the polar angle range 36° < 0^ < 144° and have the transverse energy profile 
in fhe calorimefer corresponding fo fhe profile from single elecfromagnefic shower [[TT|]. The four 
phofons are kinematically filled lo fhe Ks —)• —)• Ay hypolhesis. The dislribulion of of 

fhe kinemafic til is shown in Fig. |I] for dala and simulaled —)• KsKl evenls from fhe energy 

region ^/s = 1.04 — 1.10 GeV, in which fhe background conlribulion is small. If more lhan one 
Ks candidate is found in an even!, fhe candidate wifh fhe lowesl x^ value is chosen. For furlher 
analysis, evenls wifh x^ < 25 are selecled. 

The dominanl background process is —)■ (on^ —)• n^n^y. To suppress Ihis background 

fhe kinemafic til lo the ^ n^'nPy hypothesis is performed for events containing five or more 

phofons. Evenls with x'^o^fiy < 60 are rejected. 

The measured Ks parameters are used to calculate the mass recoiling against the Ks'. Mrec = 
^{2E - EksY - where Ek^ and pks are the Ks energy and momentum, respectively. The Mrec 
spectra for different energy intervals are shown in Fig. ^ The two peaks in the spectra for energies 


- 3 - 















Figure 3. The simulated distribution of the radius of the disappearance point (due to decay or nuclear 
interaction) for Ki mesons with E — 0.51 GeV. The cross-hatched histogram shows the contribution of Ki 
decays. The diagonally hatched areas in the zones II and IV show the contribution of nuclear interaction in 
Nal(TI) crystals. The unhatched areas in the zones I, III and V represent the contribution of nuclear interac¬ 
tion in the detector material before the calorimeter, between the calorimeter layers, and after the calorimeter, 
respectively. The interaction in the iron absorber is shown by the histogram with vertical hatching. The 
deviation of the distribution for Kl decays from the exponential near zero is due to the finite size of the 
beam interaction region along the collider beam axis (radius is measured from the geometric center of the 
detector). 


y/s > 1.06 GeV correspond to the reactions —)• KsKj^ and —)• KsKlJ. The latter reaction 

is dominated by the radiative return to the 0(1020) resonance, i.e. the process 07- In 

the following analysis we use events with the recoil mass close to Ki mass (0.40 < Mrec < 0.55 
GeV/c^), which come from the —> KsKl reaction. 

The distribution of selected events over energy intervals is given in Table |I[ In the energy 
range above the 0(1020) resonance > 1.04 GeV) 2606 events are selected. The background 
sources were studied in detailed in our earlier work ||^. The expected number of background 
events at > 1.04 GeV is 73 ±7 (36 from the beam background, 25 from nnd 12 

from —)■ (OTT^). 


4. Extraction of the Ki nuclear interaction length 

The Ki decay length Adec varies from 3.4 m at E = 0.51 GeV to 15.2 m at E = 0.1 GeV and far 
exceeds the calorimeter outer radius, about 0.7 m. The nuclear interaction length Aint is close to the 
thickness of the calorimeter sensitive volume Lcai = 34.7 cm. Therefore, a significant part of Ki 
mesons undergoes a nuclear interaction in the calorimeter. Charged particles and photons produced 
in the interaction give an energy deposition in Nal(Tl) crystals near the Ki path. Events selected 
using the criteria described above are divided into two classes, with four and with five or more 
photons. The laffer class is dominafed by evenfs, in which fhe mesons eifher decay or inferacf 
inside fhe defecfor volume. 
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Figure ^ shows the simulated radial distribution of the disappearance point (due to decay or 
nuclear interaction) for Ki mesons with E = 0.51 GeV. It is seen that the main mechanism of the Ki 
disappearance inside the detector (zones I-IV in Fig. is a nuclear interaction in Nal(Tl). From the 
experiment we extract the ratio of the number of events with five or more photons {Nsy — A^ 5 y,bkg) 
to the total number of events {N — N\±g = A ^57 + A ^47 —A^s^.bkg — ■^ 47 ,bkg )5 where A ^47 {N^y) is the 
number of selected events with four (five or more) photons, and A(t 7 ,bkg (^ 57 ,bkg) is the esfimafed 
number of background evenfs wifh four (five ore more) phofons. 

The number of evenfs wifh five or more phofons can be calculafed as follows: 

i=\ 

57 

n = wiCj 

m = >vii(l -wi)ej/, 

'"III = VVlll(l — >Vl)(l — Wll)£jj[, 

riv = >viv(l — wi)(l — >vii)(l — 

rv = (1 — wi)(l — wii)(l — wiii)(l — wiv)£v^, 

where Nq is fhe number of —)• KsKi dafa evenfs. The probabilities of fhe disappearance 

in fhe Nal(Tl) (zones II and IV in Fig. are equal to wnjv = 1 — exp(—Ln,iv/'i')> where 1/A = 
1 /Aint + l/Adec^ A// is fhe fofal fhickness of fhe l-sf and 2 -nd calorimeter layers, and Ljy is fhe 
fhickness of fhe fhird layer. The Adec is calculated as an average over defected evenfs from fhe 
energy inferval under sfudy. 

The probabilifies of decay and nuclear inferacfion before fhe calorimeter wi and befween 
fhe calorimeter layers wm, as well as fhe probabilities for an even! fo have five or more phofons 
due fo fhe disappearance in fhe /-fh defecfor zone are obfained using MC simulation. The zone 
V in Fig. ^ confains Kl decays or inferacfions oufside fhe calorimeter. Even in fhe absence of fhe 
Ki signal an even! can be selected fo fhe five-phofon class, for example, because of fhe spliffing 
elecfromagnefic shower of one of fhe phofons from Ks decay. 

The same formulae wifh fhe replacemenf of by are used fo calculate fhe number of 
four-photon evenfs. From fhe comparison of fhe predicfed rafio wifh fhe ratio 

obfained in dafa, fhe inifial value of Aint is obfained. The inifial value is fhen correcfed to fake info 
accounf fhe following effecls. 

Kl elastic scattering. The Ki elastic scaffering leads to increase of fhe effective calorimeter 
fhickness and approximately fhe same increase of fhe inferacfion lengfh. The correction fo Aint 
varies from 7.4% af E = 0.51 GeV up fo 1.9% af E = 0.69 GeV. The systematic uncertainly of Ihis 
correclion is esfimafed by variation of fhe elastic scaffering cross section used in MC simulation by 
30% and is found to be 2%. 

Beam-generated spurious photons. Abouf 10% of dafa evenfs confain addifional spurious pho¬ 
fons from fhe beam background. The spurious phofons lead to fhe Iransilion of four-phofon evenfs 
fo fhe five-phofon class and decrease fhe measured Aint value by abouf 10%. This effecl is faken info 
accounf in MC simulafion: fhe beam-background evenfs recorded during experimenf wifh a special 
random frigger are merged wifh simulated signal evenfs. The e^e^ —)■&)—)• 71^7 evenfs are used fo 
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Table 1. The energy interval (i/s), integrated luminosity (/L), number of selected events {N), number of 
background events (A^bkg), number of events with bve or more photons (Nsy, A^ 5 y,bkg), and the measured Ki 
nuclear interaction length in Nal(Tl) (Aint). The hrst error in Aint is statistical, the second systematic. 


GeV 

IL, nb ^ 

N 

^bkg 

Nsy 

^Sy.bkg 

Aint, cm 

1.02 

486 

30521 

153 

23840 

127 

31.4±0.4±1.6 

1.04 

70 

245 

3.0 

183 

2.1 

42.1±8.8±2.1 

1.05 

264 

556 

7.2 

433 

5.9 

37.4±6.1±1.9 

1.06 

432 

656 

6.6 

514 

4.7 

37.5±5.9±1.9 

1.07-1.08 

669 

516 

8.9 

391 

5.1 

45.6±7.3±2.6 

1.09-1.10 

531 

206 

5.8 

156 

4.2 

37.1±7.6±1.9 

1.11-1.13 

508 

118 

6.9 

82 

4.4 

43.3±10.3±2.2 

1.14-1.16 

674 

101 

5.9 

73 

4.1 

35.7±7.9±1.8 

1.18-1.21 

1157 

76 

7.9 

51 

5.4 

48.3±13.8±2.4 

1.22-1.29 

1916 

83 

9.8 

55 

6.7 

54.8±15.7±2.8 

1.30-1.38 

2978 

49 

10.8 

35 

7.7 

45.7±21.6±2.4 


test that the simulation reproduces well the photon multiplicity distribution observed in data. The 
systematic uncertainty in the Aint due to spurious photons is estimated to be less than 1%. 

Nonmonochromaticity of mesons. Radiative corrections, mainly due to the radiative return 
to the 0-meson resonance, lead to the deviation of the Ki average energy from ^fsjl. This effect 
is especially significant in the range =1.06 — 1.14 GeV, in which the energies of mesons 
from the reactions —)■ K^Ki^ and —)• 07 —>• KsKi^y are already sizably different, while the 
peaks from these reactions in the M^c spectra (Fig. ^,c) are not well separated. The value of the 
nuclear interaction length measured in energy interval is related to true Aint,; values for 

intervals with lower energies: 

(4.1) 

,/=l 

where F,/ is the probability for Ki produced at beam energy from interval to have (due to initial 
state radiation) energy in /* energy interval. The boundaries of the intervals listed in the first 
column of Table |] are expanded to provide full coverage of possible Ki energies, e.g. the value 
of Aint at ^/s = 1.06 GeV are used in Eq.(p~I|) for the range 1.055-1.065 GeV. The coefficients Pij 
are determined from the AA energy spectra obtained using MC simulation for different ^/s. The 
values of diagonal coefficients Pa are about 0.5 for ^/s = 1.06 — 1.14 GeV and then increase to 
about 0.9 at ^/s >1.2 GeV. The system of the linear equations ( |4.1[ ) is solved to determine Aint,;- 
The systematic uncertainty in Aint due to the nonmonochromaticity is estimated to be 1% at 
= 1.02 GeV, then increases reaching a maximum of 3% at ^/s = 1.08 GeV, and then decreases 
to 1% at 1.38 GeV. 

The procedure of Aint measurement described above is tested using simulated —)• KsKi{y) 

events. The found nuclear interaction lengths are consistent with the values Aj^^ used in simulation 
within the statistical uncertainties. The same result is obtained with A;^*^ varied by 30%. 

The corrected values of Aint obtained in data are listed in Table The first error is statistical, 
the second systematic. The sources of the systematic uncertainty are listed in Table Some of them 
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Table 2. The sources of the systematic uncertainty on Ajnt (%) for three Kl energies. 
Source 0.52 GeV 0.54 GeV 0.69 GeV 


Ki elasfic scaftering 

2.0 

2.0 

2.0 

Kl nonmonochromaficify 

1.0 

3.0 

1.0 

Spurious photons 

1.0 

1.0 

1.0 

Kl inferacfion ouf of Nal(Tl) 

3.4 

3.4 

3.4 

Background subfracfion 

1.5 

1.5 

1.5 

Defection efficiency 

2.2 

2.3 

2.6 

Tofal 

5.0 

5.8 

5.2 



E 


cal 


W7 


Figure 4. The distribution of the normalized energy deposition in the calorimeter for data events (points 
with error bars) and simulated —>■ KsKl{y) events (histogram) at y/s = 1.02 GeV. 


are described above. The uncertainty due to imperfect simulation of Ki nuclear interaction in the 
detector material before, after, and between the calorimeter layers is estimated by variation of the 
cross section used in the MC simulation by 30%. The uncertainty due to the background subtraction 
is determined by the uncertainties on the numbers of estimated background events in four- and 
five-photon classes. The systematic errors on the detection efficiencies, and ef^, associafed 
wifh Ks reconsfrucfion are cancelled in fhe rafio A^ 5 y/(A^ 5 y-|-A^ 4 y). Remaining sysfemafics is due 
fo imperfecf simulation of fhe defecfor response fo Kl mesons. The disfribufion of fhe normalized 
energy deposifion in fhe calorimeter Eca\/^/s for selecfed dafa and simulated evenfs af = 1.02 
GeV is shown Fig. Good agreemenf befween dafa and simulation specfra is seen. The narrow 
peak af £'cai/\/^ ~ 0.5 corresponds fo evenfs wifh undefecfed Ki meson. Since fhe fhreshold on fhe 
Kl energy deposifion is low (20 MeV), we do nof expecf any significanl systematic uncerfainfy due 
fo inaccuracy in fhe simulation of fhe defector response. The uncerfainfy in Table ^ associafed wifh 
fhe defection efficiency is determined by simulation sfafisfics. 
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Figure 5. The nuclear interaction length in Nal(Tl) as a function of the Kl energy. Points with error 
bars represent the result of this work. The solid curve shows the nuclear interaction length used in the 
UNIMOD [P simulation package. The dotted line is the same dependence used in the GEANT4 package, 
version 9.5 (low-energy physics) [^. 


5. Conclusion 


Using kaons from the reaction —)• KsKi we have measured the nuclear interaction length 

in Nal(Tl) in the Kl energy (momentum) range from 0.51 (0.11) to 0.69 (0.48) GeV (GeV/c). The 
data were collected by the SND detector at the VEPP-2M collider. We did not find any other 
measurements of the Ki^ nuclear length at such low momenta in literature. 

The energy dependence of the Ki nuclear interaction length measured in this work is shown in 
Fig. ^ in comparison with the nuclear lengths used in the UNIMOD ||I|] and GEANT4 packages. 
The nuclear length for GEANT4 is obtained using test simulation, in which Ki mesons interact in 
a large block of Nal(Tl). The is calculated from the probability of disappearance in a thin 
layer of material. 

Our data are in good agreement with UNIMOD. The values of the nuclear length used in 


GEANT4 version 9.5 (the model of hadron physics FTFP_BERT [12]) contradict to the measured 


values at energies below 0.55 GeV. The results of this work may be used to refine the A^ nuclear 
cross section in GEANT4. 


Acknowledgments 

This work is supported by Russian Science Foundation (project No. 14-50-00080). 

References 

[1] A. D. Bukin et ai. Preprint INP 86-18, Novosibirsk, 1986. 

[2] R. Brun, R Bruyant, M. Maire, A. C. McPherson and P. Zanarini, CERN Report No. 
CERN-DD/EE/84-1, 1987. 

[3] S. Agostinelli et al, Nucl. Instrum. Methods Phys. Res., Sect. A 506, 250 (2003); J. Allison et al, 
IEEE Trans, on Nucl. Science 53, 270 (2006). 

- 8 - 




















[4] A. N. Skrinsky, “VEPP-2M status and prospects, and 0-factory project at Novosibirsk,” in Proc. of 
Workshop on physics and detectors for DA<t>NE 95, 1995, vol.IV, p.3. 

[5] M. N. Achasov et ai, Nucl. Instrum. Methods Phys. Res., Sect. A449, 125 (2000). 

[6] M. N. Achasov et al, JETP 103, 720 (2006). 

[7] E. A. Kuraev and V. S. Eadin, Sov. J. Nucl. Phys. 41, 466 (1985) [Yad. Eiz. 41, 733 (1985)]. 

[8] G. Bonneau and E. Martin, Nucl. Phys. B 27, 381 (1971). 

[9] A. M. Makhov, Preprint INP 92-66, Novosibirsk, 1992. 

[10] K. Haenssgen et al. Preprint KMU-HEP, 80-07, Leipzig, 1980. 

[11] A. V. Bozhenok, V. N. Ivanchenko and Z. K. Silagadze, Nucl. Instrum. Methods Phys. Res., Sect. 
A379, 507 (1996). 

[12] Geant4 Physics Reference Manual, in http://cern.ch/geant4/support/userdocuments.shtml [cern.ch] 


- 9 - 



